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Abstract 

We consider in this paper a spray constituted of an incompressible viscous gas and of small droplets 
which can breakup. This spray is modeled by the coupling (through a drag force term) of the incom- 
pressible Navier-Stokes equation and of the Vlasov-Boltzmann equation, together with a fragmentation 
kernel. We first show at the formal level that if the droplets are very small after the breakup, then the 
solutions of this system converge towards the solution of a simplified system in which the small droplets 
produced by the breakup are treated as part of the fluid. Then, existence of global weak solutions for 
this last system is shown to hold, thanks to the use of the DiPerna-Lions theory for singular transport 
equations. 

1 Introduction 

Sprays are complex flows which are constituted of an underlying gas in which a population of droplets (or 
dust specks) are dispersed, cf. |18| . There are various possibilities for modeling such flows, depending in 
particular on the volume fraction of the liquid phase (cf. [H] for example). 

We focus here on the case when the volume fraction occupied by the droplets is small enough to be 
neglected in the equations (such sprays are called thin sprays, cf. [H]), so that the modeling of the liquid 
phase can be performed by the use of a pdf (particles density function) which solves Vlasov-Boltzmann 
equation (cf. [Hid]). Denoting / := f(t,x,£,r) > the number density of droplets of radius r which at 
time t and point x have velocity £, the Vlasov equation writes 

&/ + €-Vx/ + v € -(/r) = Q(/), (l) 

where T represents the acceleration felt by the droplets (resulting from the drag force exerted by the gas), 
and Q is an operator taking into account the complex phenomena happening at the level of the droplets 
(collisions, coalescences, breakup). 

We also restrict ourselves to the case when the gas is incompressible and viscous, which is for instance 
the usual framework when studying the transport of sprays in the upper airways of the human lungs, cf. 
[T7] . Accordingly, the gas is modeled by the incompressible Navier-Stokes equation 

V x ■ u = 0, (2) 
Pg [d t u + V x • (u $ u)] + V x p - pA x u = F ret , (3) 

where p g is the constant density of the gas, u := u(t,x) 6 M 3 is its velocity, p is its constant (dynamic) 
viscosity, and F ret is the retroaction of the drag force: 
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F ret (i,x) = - / / ^p ir 3 fTd£dr. (4) 

JO JR 3 6 

Finally, pi is the constant density of the liquid (so that the mass of the droplets of radius r is | pi r 3 ). 
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For particles with small Reynolds numbers, the drag force is given by the simple formula (known as 
"Stokes' law") 

r(*,x,£,r) = -- , 5 

2 pi r z 

that we shall systematically use in the sequel. 

The modeling of the breakup phenomena is important in the applications and has led to various models 
appearing in the literature (cf. [T]). A typical form of the breakup kernel is obtained when assuming that 
the droplets after breakup have the same velocities as before breakup. The operator writes then 



Q(/)(t,x,€,r) 



b(r*,r)v(£,r*)f(t,x,$,r*)dr*, 



(6) 



where v := r) > is the fragmentation rate, and b := b(r* , r) > is related to the probability of ending 
up with droplets of radius r out of the breakup of droplets of radius r* . We finally obtain the following 
system 

V x • u = 0, (7) 
S t / + C-V x / + V € -[/r] = Q(/), (8) 



Pg 



d t u + V x • (u ® u) + V x p - /iA x u 



(9) 



Because of the dependence w.r.t. r of T in cq. ((SJ), we see that the drag force acting on the small droplets 
leads to very quick equilibration of their velocity with the velocity of the gas. It is therefore natural to try to 
write a set of equations replacing (jl]) - © , in which the small droplets are considered as part of the gaseous 
phase. 

Indeed, in the context of the numerical simulation of cq. |T]) thanks to a particle method, the small 
droplets which are produced because of the breakup can lead to a high computational cost. Once integrated 
(at the level of the continuous equations) in the gaseous phase, they can be discretized along with the gas 
thanks to a finite volume scheme, and their contribution to the computational cost remains in this way 
reasonable. 

In order to perform a mathematical study of the system obtained by such an approximation, we consider 
the simplest possible case, namely when the fragmentation rate v takes the form v{£,r) := r l r > r2 , for 
some constants r 2 > and r > 0, the latter being seen as a characteristic time of fragmentation. We also 
assume that the aerosol is bidispcrsed: only two possible radii r\ > T2 exist for the droplets, and the result 
of the breakup of particles of radius r± are particles of radius T2 ■ This implies that the density / splits in 
the following way 

/(*,*)£>>") = /l(*,X,£)5 r =ri +/2(£,X,£)£r=r 2! 

and that system ((3J) - © rewrites, after having normalized all the constants (except ra): 



V x • u = 0, 

dtfi + € • Vx/i + V c ■ [fx (u - £)] 
dth + i ■ V x / 2 + V C 



h 

„3 ' 



8 t u + V x • (u ® u) + V xP - A x u = - / (u - £) h d£ - r 2 / (u - £) f 2 d£. 

JR3 JS.3 



(10) 

(11) 

(12) 
(13) 



As we already explained, we are interested in the asymptotic regime when the particles resulting from 
breakup are becoming smaller and smaller. In our bidispcrsed model, this reduces to study the limit r 2 — > 0. 
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Denoting p = / fi d£, and noticing that 



el 

dt 



/^d^dx + r 3 , / / / 2 M_d€dx- 



dx 



T 3 



/ f /i|^-u| 2 d€dx 

Jt 3 7 k 3 



+ 



r 3 / 2 |£-u| 2 d£dx 



|V x iiMdx = 0, 



r 2 v/t 3 Jr 3 Jt 3 
we see that (at the formal level) rf /2(i, x, £) — > p(t, x) ^ =u (t )X ) when r 2 — > 0. 
Integrating cq. (fT2"|) against r| d£, we end up with 

d t p + V x - [pu] = f hdt (14) 

Then, integrating eq. (fT2"|) against r| £d£ and adding the result with eq. (fT"3"|) . we obtain 

d t ((l + p)u)+V x -((l+p)u®u)+V x p-A x u = - / ( u -0/id£+ / Mde 

JR 3 JR 3 

Combining this last equation with eq. (|14[) and replacing the notation /i by /, we write down the system that 
we wish to study (we close it with periodic boundary conditions for the mathematical proof of existence): 



V x u = 0, 

%> + V x -[pu] = / /d£, 

,/R 3 

9e/ + €-V 3t / + V*-[(u-£) /] = -/, 

(l + p)[<9 t u + V x -(u®u)l +V x p-A x ii = 2 / (£-u)/d£, 



(15) 
(16) 
(17) 
(18) 



where p := /?(t,x) > 0, u := u(*,x) G M 3 , p := p(i,x) > 0, / := /(t,x,£) > 0, and i > 0, x G T 3 , £ G P 3 . 
Let us recall that p represents here the "added density" resulting from the very small particles. This is the 
reason why, though we normalized the fluid density in (|13p by p g = 1, we have in (|18p the term (1 + p) in 
front of the convective part of the fluid equation. For a more detailed version of the previous computation, 
see [S]. 

Our goal is to study the existence theory for this system, completed with the following initial data: 

/(0,x,O = /in(x,O>0, xeT 3 , £eR 3 , (19) 
p(0,x) =pta(x) >0, xeT 3 , (20) 
u(0, x) = U;„(x) G T 3 , xeT 3 . (21) 

Before stating our result, let us introduce a few notations. If there is no ambiguity on the considered time 
interval, we simply denote by Lf(L«) and L?(L«^) the spaces L p ([0, T]; L 9 (T 3 )) and I/([0, T]; L 9 (T 3 x R 3 )), 
for any pair of exponents (p, q) £ [1, oo] 2 . In particular when p = q we simply use the notation L p . We adopt 
the same convention for Sobolev spaces Wi' p and H™ (m € N) and denote by H" m (or equivalently R~ m (T 3 )) 
the dual of the latter. When the subscript " div " is added to any space, the corresponding subspace are 
the divergence-free (in x only) elements of the ambient space. We denote by P the Leray projector onto 
divergence-free vector fields 

P = Id - V X A" 1 div x , 
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easily defined on @ , ([0,T] x T 3 ) thanks to Fourier analysis, extended by duality to ®'([0, T] x T 3 ). 

Finally, if h is a scalar function defined on M. + x T 3 x R 3 , we define the following moments for h := h(t, x, £) 
(and for a£t) 



m a h(t, x) 



/ |£| a Ji(t,x,e)d& M a h(t):= [ m a (/i)(i,x)dx, miA(t,x) := / d t 

JR3 ,/T 3 JR 3 



Thanks to those notations, we are able to write down our main Theorem: 
Theorem 1.1. Let T > and assume that 

/to G L°°, (1 + |£| 2 )/ in G L\p in G L°°,Uto G I4 V . (22) 

TTien £/ie system (|15|) - (|2ip admits a global weak solution (p > 0, / > 0, u G R 3 ) smc/i i/iai 

p G L°°([0,T],L 5/3 (T 3 )),/ G L°°([0,T] x T 3 x R 3 ),u G L 2 ([0, T], H^ iv (T 3 )) n L°°([0, T], L 2 (T 3 )). (23) 
Moreover, the triplet (p, /, u) satisfies the following energy estimate 

UM 2 f(t) + \\^/TTWju(t)\\l 2{T3) }+ r||Vxu( S )|| 2 2(T3) d S + | f / |u - £| 2 /d£dxd s 

< i {M 2 / in + HVl + Pin u in || 2 2(T3) }, (24) 

and £/ie bound 

H/||l°°([0,T]xT 3 xR 3 ) < e 2T ||/in||L°°(T 3 xR 3 )- 

Remark 1.1. We explain here the meaning of " weak solutions " in the above Theorem. For p G [l,oo], 
when u G L p , / G L°° and p G L p , equations (jl5|) - (|16|) /lave a cZear meaning in the distributional sense. 
However the meaning of (|18|) is noi completely obvious, and we rewrite this equation, using both (|18p and 

OH) , as 

p{<9,[(1 +p)u] +div x (l + p)u®u } - A x u = P{2mi/-um /}. (25) 

For £/ie initial condition, we adopt the following general definition : consider v in G ^'(T 3 ) and id in 

Lj oc ([0,T];H _m (T 3 )) for some integer m G N. A distribution v G - oo,T[xT 3 ) is a solution of the 
Cauchy problem 

d t v = w, v (0)=v ln , 
if the support of v is included in R + and if 

d t v = w + 5 Q ®v- m in 3>' (} -oo,T[xT 3 ), 

where w is the extension of w by on R_ . 

The study of the existence of solutions to coupled (through drag force interaction) fluid-kinetic equations 
is now a well-established subject. 

As far as viscous equations are concerned for the modeling of the gaseous phase, we would like to quote 
the works on the Vlasov-Stokcs equations in [T^], on the Vlasov/incompressible Navier-Stokes equations 
in [2| El [50], on the Vlasov-Fokker-Planck/incompressible Euler equations in [5] , on the Vlasov-Fokker- 
Planck/incompressible Navier-Stokes equations in [T5], and on the Vlasov-Fokker-Planck/compressible Navier- 
Stokes equations in [IT)] . 

Without any viscosity or Fokker-Planck damping, the study is more difficult, and only smooth local 
solutions are known to exist, in the compressible setting (cf. [3] and [!%])■ 
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Our system (fT5j) - (fT8"|) is reminiscent of the Vlasov/incompressible Navier-Stokes equations with a vari- 
able density. It is known that variable densities in the incompressible Navier-Stokes equations lead to extra 
difficulties (w.r.t. constant densities), cf. (7J[I3]. Here extra difficulties (w.r.t. [5]) appear in the final 
passage to the limit (that is, in the stability result for solutions satisfying the natural a priori estimates of 
the problem). They are linked with nonlinearitics which are specific of our model, and which necessitate a 
careful treatment using refined versions of Lemmas presented in [10] . 

The equations coming out of the theory of sprays are known to be difficult to approximate in a good way 
(that is, in a way in which the a priori estimates satisfied by the equations themselves, such as the energy 
estimate, are also satisfied or well approximated by the approximating equations). This difficulty appears 
in this paper, and the approximating scheme which is used is rather complicated and necessitates successive 
steps, which are reminiscent of those used in [5J. 

Our paper is structured as follows: Section 2 is devoted to the Proof of existence of solutions to an 
approximated version of this system. The approximation is removed in Section 3 and Theorem ll.il is proven 
there. Finally, we present in a short Appendix some autonomous results which are used in the Proof of 
Theorem 11.11 

2 Existence for a regularized system 

In this section, we begin the Proof of Theorem 1 1.1 1 bv introducing a regularized system. 
2.1 Definition of the regularized system 

In order to do so, we introduce a mollifier ip £ <G ^""(T 3 ) for all e > 0. We also introduce a truncation in 
this way: we define 7 e a nonnegative element of S?(R 3 ), whose support lies in B(0, 2/e), bounded by 1, and 
equal to 1 on the ball B(0, 1/e). 

Our regularized system writes 



with initial conditions f? n G ^(T 3 x K 3 ), pf n G < rf oc (T 3 ) and uf n G ^(T 3 ) approximating respectively f in 
in all Lf oc (T 3 x R 3 ) (p < oo) and in L^T 3 x M 3 , (1 + |£| 2 )dxd|), Pin in all L P (T 3 ) (p < oo), and u in in 
L 2 (T 3 ). 



Next subsection is devoted to the Proof of existence of solutions to this approximated system (that is, 
for a given parameter e > 0). 

2.2 Existence of solutions for the regularized system 

In this Subsection, the parameter e is fixed, and we drop out the corresponding indices in order to make the 
formulas more readable. 

We fix a triplet of initial data (ui, h, v) £ x ( < ^ 1 n L°°) x < ^ 71 (the two last ones being nonnegative). 



div x u e = 0, 

d t pe + (U s * <p e ) ■ V x Pe = mo(/e7e), 
( dtfe + $, ■ V x /e + V € ■ [(U E * <p e - £)/=] = -A, 

P{(1 + Pe)[d t u 6 + (u e *<p e ) ■ V x u e ] } - A x u £ = 2P{m 1 (/ e7e ) 



— u e 



(26) 
(27) 
(28) 

(29) 



5 



The considered system (in which we do not explicitly write the initial data) is then the following: 

'div x u = 0, (30) 
t?t/9+(u*p)-Vx/0 = mo(/7), (31) 
dtf + i- V x / + V € • [(u * cp - £)/] = -/, (32) 

P{(1 + p) [d t u + (u*v>) • V x u] } - A x u = 2P{ mi (/ 7 ) - um (/ 7 )}. (33) 
We are able to prove the 

Proposition 2.1. For any triplet of initial data (u>, /i, i/) £ x ( < ^ 51 n L°°) x "if 1 (the two last ones being 
nonnegative) , there exists a weak solution (in the sense of Remark ] 1 . 1}) to the system h3(M 



Proof of Proposition [2TD 

We shall show existence of a solution to the nonlinear system ([3L)p - (|33|) by applying Schauder's fixed 
point Theorem to the following mapping: 

S : ^(^div) ► ^?( L div) 

u^u°, 

where u° will be built thanks to the following steps : 

(i) We first consider /° £ ff 1 , unique classical solution of 

d t f + 1 • v x r + v € • [( U * tp - on = -r, m 

with h as initial datum, 

(ii) then p° £ c € x is defined as the unique classical solution of 

d tP « + (u* <p) ■ V x p° = m (r 7 ), (35) 

with v as initial datum, 

(iii) and finally u° £ ^ t (Lj iv ) is built as the unique divergence- free weak solution of 

P j (1 + p«) [d t u« + (u* <p) ■ V x u«] | - A x u* = 2P{ mi (/*7) - u*m o (/ 7)}, (36) 
with ui as initial datum. 

In the sequel Pt(^, h, v) will denote any nonnegative function (which also may depend on e, but as explained 
before we omit the corresponding index here) that splits into a finite sum of positively homogeneous functions 
of strictly positive degree of one of the following terms : HwHh 1 , Halloo or |H|oo- Pt(w, h,v) may change 
from one line to another, but will always have the structure that we just described. 

In the next paragraph, we show that /°, u° and p° are well defined. 
2.2.1 Existence and uniqueness of /°, p° and u* 

We know that u*ip £ ^(^€^) and V x (u*<^) £ L°°, so that the characteristic curves of and ([33]) arc 
globally well-defined. We hence easily obtain the existence and uniqueness of two nonnegative functions /° 
and p° (classically) solving (|34|) and (|3"5|) . Thanks to the maximum principle, we also know that 

uriioo < e 2T ||^IU, 
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from which we get 

||mo(r7)lloo + ||m 1 (r 7 )|| 00 <C 7 e 2T ||/ i || 00 <P T (a;,ft^), (37) 

and hence 

||p°||oo<Pt(w,M). (38) 

We now consider /° and p° as given. The existence and uniqueness of u° E H| x C ^(L 2 ), divergence- 
free weak solution of (|36|) (with u for initial data) may be obtained thanks to the usual Galerkin method. 
Since this type of construction is rather standard, we won't detail it here, and refer to [5] for a precise 
treatment of this procedure in our system, or for instance to [3J, Chapter 3, for a more generic approach. 

In next subsection, we obtain natural estimates for the quantity u°. 
2.2.2 Estimates for u° 

Taking u° as test function in (|36|) and using (|35p . we obtain the following estimate (since u*ip is divergence 
free and /° is nonnegative) 



\\j (l + /^,x))|u^,x)| 2 dx)+ f [ |V x u*( S ,x)| 2 dxd S 
z Ui 3 J Jo Jt 3 

u*( S ,x)| 2 m (/ <> 7)(s,x) dxds 



< - 






2 


+ 2 


/; 






1 




+ 2 


{ 



3 '* 



JT 



(l + Kx))Mx)| 2 dx}>, (39) 

T 3 



so that using p° > and Gronwall's lemma together with estimates (|3"7| - (f3"5)) , we get 

||u*||Lf (14) + H u lL?(Hi) < Pr(w, h, v) exp(P T (u;, h, v)). (40) 
Taking then <9tU° as a test function in ([55)) . we get (using again p° > and /° > 0) 

ll^ll 2 + \ ||V x u°(t)|| 2 < ||(1 + p o)(u *^)|| 00 ||V x uO|| 2 ||d t ul 2 + 2 !| mi (r 7 )|| 2 ||a t u°|| 2 
+ 2 ||m (r 7 )|U K|| 2 ||^u*|| 2 + \ ||V x w!| 2 , 
so that using Young's inequality, estimate (|40p above, and estimates (|37 p -(|38 p . we end up with 

||9tu°|| 2 < Pt(w,M exp(Pr(w,M)) [||u|| Lr (i4) + l] , (41) 

thanks to the elementary convolution inequality ||u*y>||oo < II U H (t-E) II ^11 2 ■ 

Given any triplet (w, /i, z/) € x "if 1 n L°° x of initial data, we can then define the mapping 

S : (Ljiv) ► ^"(Ldiv) 
U I — ► u°. 

In next paragraph, we use Schaudcr's fixed point Theorem for this mapping. 
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2.2.3 Application of Schauder's Theorem 

According to the estimates of the previous paragraph, we see that for any u <G ^f t °(Ljj t ), 

||S(u)||i 4 o (I( , ) <P r (w,/i J i/)eicp(Pr(w,ft,i/)) := R- 

Hence, S sends (the closed ball of radius R of the normed space ^ (L^ iv )) to itself. 

Thanks to estimates (|40p and (|41[) . a bounded subset of ^ t (L^ iv ) is sent by S to a bounded subset of 
Hj x , and Lemma FOl of the Appendix allows us to conclude that S(B^) is relatively compact in ^ t (Lj iv ). 

It remains to study the continuity of S. Notice first that for any u <E "^(Ljj. ), any weak bounded solution 
of or ([55)) is a renormalized solution in the DiPerna-Lions sense introduced in [TO]: the absorption 
coefficient (= 1) of the first equation and the right-hand side of second one are bounded and one easily 
checks that the vector fields (t, x, £) H> (£, u * (pit, x) — £) and (t, x) H> u*y;(t,x) satisfy all the desired 
assumptions given in [TU]. Now if (u„)„ converges to u in ^°(L^ iv ), we know from the previous step that 
(S(u„))„ e N has a converging subsequence. Showing that the whole sequence converges to S(u) reduces hence 
to prove that it has only one accumulation point, namely S(u). Assume therefore that v G ^ t °(Lj. ) is 
an accumulation point of (S(u„))„ s n and denote by a the corresponding extraction. Since (u a ^) n is still 
bounded, using (|40p - (|4ip and adding a subsequence if necessary, we may assume that v is also the limit of 
(S(u„))„ s n in Hj x , for the weak topology. On the other hand, since (u CT („)) n still converges to u in ^(L 2 ^), 
DiPerna-Lions stability result of [10] ensures that (/ (T ( ri ))n and (p a ( n ))n both converge in L^ oc (for allp < oo) 
to / and p, the corresponding (unique) weak and bounded solutions of the associated equations (with initial 
data h and v and vector field defined by u). At this point, we have enough (strong) convergences to ensure 
that all the nonlinear terms of the fluid equation indeed converge to the expected limit : v G x is hence 
a weak solution of (|36p with w as initial datum. We already mentioned the uniqueness of solutions for this 
equation, so that we eventually get v = u° = S(u). 

S is hence a continuous map from a closed convex nonempty set of a normed space to itself, and it has 
a relatively compact range. Thanks to Schauder's fixed point Theorem (see 1 1 1 j for instance), S has a fixed 
point which is a solution to (f3"0)) - (f3"3")) with initial data (w, h, v). This concludes the Proof of Proposition 

o □ 

Reintroducing the parameter e > 0, we have obtained the existence (for any T > 0, e > 0) of a divergence- 
free vector field u e <G Hj x , and two nonnegative functions f e £ ^ x ^ and p e £ solutions of system (|26J) 
- ([Ml). 

Note that since u e ★ ip E e L^. and f? n is compactly supported, one can show without difficulty that 
f E (t,-,-) remains compactly supported for a.c. t (with a support depending on e). 

In next section, we pass to the limit when e — ¥ in the functions u e , f e and p £ , and show that their limit 
is a (weak) solution of (fT5")) - (fT5|) . 

3 Passage to the limit and Proof of the main Theorem 

We keep on proving Theorem 1 1.1 1 in this Section. We start with a local result of existence: 

Proposition 3.1. Under the assumptions on the initial data of Theorem ] 1 . 11 there exists an interval J m := 
[0,t+] on which the system (|15p - (|18p admits a weak solution (p, /, u) (cf. Remark \l.l\ for a definition of 
weak solutions) . Furthermore p and f are nonnegative, and 

peL°°(J in ;L 5 / 3 (T 3 )), 

/eL°°(J in xT 3 xlR 3 ), 

u e L 2 (J in ; H^TT 3 )) n L°°(J in ; L 2 (T 3 )). 
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Proof of Proposition [37T1 : 

We begin with the Proof of bounds for the solution of - (|!?9")) . which do not depend on e. 

In the sequel P^fu^, p m , f m ) [or more simply P™] will denote a polynomial function (with positive 
coefficients) of the quantity \\y/l + p m Ui n || 2 + ||/in||oo + -^/in- Note that will always be independent of 
e, but may depend on T. 



3.1 Uniform bounds with respect to e 

Let us first recall a classical Lemma linking velocity moments, the proof of which may be found in [B] 
(Lemma 1, Section 3.3): 

Lemma 3.1. Let 7 > and h be a nonnegative element o/L°°([0, T] x T 3 x R 3 ), such that m 7 h(t, x) < +00 
for a.e. (t, x). The following estimate holds for < a < 7: 

m a h(t,x) < (J^Tr\\h(t,x,-)\\ L ^ {VL p + l^j m 7 /i(i,x)W. 

Thanks to the maximum principle in eq. (|28[) . we observe that 

IIMloo < ||/in||ooe 2T , (42) 
from which, thanks to Lemma |3. 11 we deduce 

IK/eWI! 5 / 3 <PT n A^/ £ (*) 3/5 , (43) 

llmiAWII^P^AW 4 / 5 . (44) 

Since f e is compactly supported and solves eq. (|28[) strongly, one gets by multiplying this equation by |£| 2 
and integrating in x, 

d /* 

-M 2 f e + 3M 2 / e = 2 (u £ * tp £ ) • mi(/ B ) dx. (45) 
at J T 3 

Using the previous estimates, we hence have 

=1 

^M 2 / £ + 3M 2 / £ < ||u e (t)|| 5 || mi / e (t)|| 5/4 
<P£ ||u £ (t)|| 5 M 2 / e (t) 4 / 5 , 

from which we easily deduce 

M 2 fe(t) < {m^ + P^ ||Ue( S )|| 5 d S J . 

Using the Sobolev injection H 1 (T 3 ) L 5 (T 3 ), one gets 

M 2 / e (t)<P£ jl + ^ ||u £ ( S )|| 2 ds + ^ ||V x U £ ( S )|| 2 ds| 

Using the previous estimate with (|4"3"|) - (|4"4")h we obtain 

||mo/ s (t)|| B/3 < P£ {l + j K(s)|| 2 d fl + j ||V x u £ ( S )|| 2 <L 
\\mif £ (t)\\ 5 /4 < Pr {l + j l|u E (s)|| 2 d S + j ||V x u £ ( S )|| 2 d 



(46) 



4 

s\ . (48) 
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The energy estimate ([55)1 is satisfied with v := pf n , oj := uf n , and rewrites 

\{ f (l + p £ (i,x))|u £ (i,x)| 2 dxl + f f |V x u £ ( S ,x)| 2 dxds 
1 Ut 3 J Jo Jt 3 



< 



-7; [ ( |u 6 (s,x)| 2 TO (/ e 7)(s,x)dxds 
* Jo Jt 3 

2/ / u e (s,x) • mi(/ 6 7)(s,x)dxds 
Jo Jt 3 

(l + /4(x))|uf n (x)| 2 dxl, 



so that (since p e is nonnegative) we obtain 

^IM*)lla+ j[ !|V x u £ (.s)|| 2 d.s< i|| v / T+PLuf n || 2 + 2^ Hue-mi/eWHida 

<P£ + 2 / ||u e -m 1 / e (s)|| 1 ds. 
Jo 

Using Holder and Young inequalities, we get 

^K(*)lll + \f II V x u e (.s)|| 2 ds < P*f + P£ J* ||u £ (.s)|| 2 d S + P* jT* || mi / £ ( S )|| 2 /4 d*. 
We then obtain thanks to 

KWII2 + ^ VxUe^H 2 ds < P* + P* J* |||u £ (,s)|| 2 + Vu e (a)||^ dajV 

Denote by z the maximal solution of the Cauchy problem 

z'{t) = ¥ i «z{t)\ z{Q) = ?t, 

defined on some maximal interval Jj n . Then if J; n is the closure of Ii n /2, using the nonlinear Gronwall 
Lemma |4"TTI we get for t G Ji n : 

K(t)|| 2 + /'j|Vu £ ( s )j| 2 d.s<P|| L ~ := . 
Jo 

Remark 3.1. Notice that thanks to Lemma\4-S\ (changing P™ if necessary) we always have \ Ji n \ > 1/Py- 



We deduce from the previous local estimate that (u £ ) £ is bounded in L°°(j in ; L 2 (T 3 )) n L 2 (j in ; H^T 3 )). 
Using we hence get the boundedness of (M 2 / £ ) e in L°°(J in ) and then, with (Ti5|) and (jUJ), we see that 
(m / £ ) £ and (mi/ e ) e are respectively bounded in L°° ( J in ; L 5/3 (T 3 )) and L°° ( J in ; L 5/4 (T 3 )) . Thanks to a 
classical transport estimate, using (|28p. we see that (p £ ) £ is bounded in L°° ( J in ; L 5/3 (T 3 )) . From (|46p . we 
deduce that for all e > and t G Jj n , 

M a / e (t) < (49) 

3.2 Compactness properties 

Recalling Proposition 12.11 of existence of a solution to the regularized problem (for a given e), we deduce 
from the previous (uniform in e) bounds the existence of p £ L°°(Ji n ; L 5 ^ 3 (T 3 )), / <G L°°( Jj n ; L°°(T 3 x 
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such that, up to a subsequence, 

( Pe ) £ - •oinL~(J ln ;L B / 3 (T 3 ))-*, 

£— >0 

(/ e ) E - / in L^J^L 00 ^ 3 x W 3 )) - *. 

Using the boundedness of (A/ 2 /e)e in L°°(Ji n ) and adding a subsequence if necessary, one manages to also 
show that 

(m / s ) e and (m (/ e7e )) £ - m / in L°°( J in ; L 5 / 3 (T 3 )) - *, 

e— >0 

(mi/,) t and (mi(/ e7e )) e - mi/ in L°°(J in ; L 5 / 4 (T 3 )) - *. 

Finally, using the bounds on (u e ) E and the Sobolev injection H 1 (T 3 ) L 6 (T 3 ), we get the existence of 
u G L°°(J in ;L 2 (T 3 )) nL 2 (J in ;Hi iv (T 3 )) such that 



(u e ) E ^uinL^J^H^T 3 )), 
(Ue) 8 - uinL 2 (j in ;L 6 (T 3 )), 
(u £ ) £ ^ o U in L°°(j in ;L 2 (T 3 )) - *, 

and of course (u £ *<p £ ) £ converges also towards u for the same topologies. 
3.3 Weak convergence of nonlinear terms 

In view of the previous weak convergences, it only remains to check that the nonlinear terms converge to the 
expected limits, so that the weak limit triplet (/, p, u) will indeed be a solution of our system. Notice that 
since all the bounds and weak extractions are performed on the local interval Jj n , the constructed solution 
will only be local in time. In the last subsection, we shall explain how to extend it. Until then, through all 
the current subsection, the index t will refer to the intervall Ji n in the notation "Lf (E)". 

Notice that the system can be written 

div x u £ = 0, (50) 

d t /e + div x4 (a £ / £ )-2/ £ = 0, (51) 

d t pe + div x (p £ (u £ *tp e )) = mo(/c7e)j ( 52 ) 

P{$[(1 + p £ )u £ ] + div x [(l + p £ )(u £ * <p e ) ® u £ ] } - A x u £ = p{2rm(/ £7£ ) - u £ m (/ £7£ )}, (53) 

where a £ (i,x,|) := (£, [u E (t) * (f e ](x)). 

In order to handle the nonlinear terms, we shall use several times the Proposition 14. II of the Appendix. 
We denote by ^ s the vector space of bounded measures on J m x T 3 . 

All the coming facts and their proofs are true up to some (finite number of) extractions that we don't mention 
in the sequel. 

Fact 1 : The products (p e (u £ * tp £ ))e an d (pe u e)e both tend to pu in ^0 S . 



Proof of Fact 1: Using Proposition 14.11 of the Appendix, these two terms are handled in the same 
way: the velocity term (convoluted by ip e or not) plays the role of (a £ ) £ , and is bounded in L 2 (W X 2 ), 
whereas (p e ) e plays the role of (6 £ ) £ and is bounded in L t °°(L x /3 ) ((5/3)' = 5/2 < 2* = 6). Thanks to 
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(f52"| and the previous bounds (see subsection 13. (dtp £ )e is bounded in L 2 (H x m ), (where m is taken 
large enough). We hence have 

(p e (u e * ip £ )) £ pu, in M s - w*, 
(/5 e u e ) e — pu, in — w * . 

Fact 2: ((p e + l)u e , u £ ) 2 — > ((p + l)u, u) 2 . 
Proof of Fact 2: We wish to prove that 

lim / / (p c - + l)|u e | 2 dxdi = / / (p+ 1) |u| 2 dxdt. 

Jj in JT 3 J J in JT 3 

First write, since u e is divergence-free, 

/ / (p E + 1) |u e | 2 dxdt = f [ (p E + l)u £ • u £ dxdt 

J J in JT 3 ij in JT 3 



(54) 



/ / P[(ft + 1K] -u £ dxdt. 

J J in JT 3 



Harmonic analysis and singular integral theory allow to show that P is bounded from (L x ) to itself, 
for all p €]l,oo[ and q £ [l,oo] (see for instance |19|). Since the strong continuity of an operator 
implies its weak sequential continuity, we see that P is sequentially continuous from L 2 (L 30 / 23 ) to 
itself, equipped with the weak topology. But we have seen in Subsection 13.11 that (u £ ) e is bounded 
in L 2 (H X ), which is embedded in L 2 (L X ) by Sobolev injection, and {p e ) e is bounded in L^°(L^ 3 ) so 
that by Holder inequality, ((p £ + l)u e ) e is bounded in L 2 (L X ^ 23 ) and hence (up to a subsequence) 
(F[{pe + 1) Uc] J converges weakly to P[(p + 1) u] in L t 2 (L 30/23 ). We now can use Proposition K2~\ of 
the Appendix, with q = 2 and r = oo. Indeed, 

- (u e ) e is bounded in L 2 (W X ' 2 ) n L t °°(L 2 ), 

- since d t P = Fd t , J53J implies that (d t F[(p E + l)u e ]) is bounded in ~L\ (H x m ) for m large enough, 

- (30/23)' = 30/7 < 2* = 6. 

Hence the product |p[(p e + 1) u E ] • u e | converges to P[(p + 1) u] • u in M s — *• In particular, using 
lji„xT 3 as a test function, we get Fact 2. 

Fact 3 (u e ) E and (vL e -kip e ) e both strongly converge to u in L| x . 
Proof of Fact 3: We estimate 

/ / |u s -u| 2 dxdt< / / (1 + p £ ) |u £ - u| 2 dxdi 

J J in JT 3 J Jin "'T 3 

= / / (l + ( o e )|u e | 2 dxdt + / / (l + p e )|u| 2 dxdt 

J J in JT 3 J J in JT 3 

-2 1 / (1 +p e )u e ■ udxeh. 

J Ji„ JT 3 

The first term of the second line converges to the same expression, but without e: this is exactly 
Fact 2 proven above. We have the same behavior for the second term of this line: (p £ ) £ converges 
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weakly to p in L^°(L^/ 3 ) — * and |u| 2 G (L x ). Eventually, for the third and last term, we use the 
(already proven) convergence of [(1 + p £ )u £ ] £ towards (1 + p) u in L 2 (L X ^ 23 ), and the embedding 
u G L 2 (L x ) L 2 (L x °^ 7 ). Strong convergence of (u E ) E is then directly transfered to (u E * ip e )e- 

Fact 4 (u E ) E converges strongly to u in all Lj(L x ), /or c < 2 and d < 6. 

Proof of Fact 4: Thanks to Fact 3, we have strong convergence of (u e ) e , and hence almost everywhere 
convergence. Since this family is bounded L 2 (L x ), we get the desired convergences. 

Fact 5 (/ E (u E * tp s )) e and (u E mo(f £ j £ )) £ converges weakly in L 1 to respectively f u and umj(/). 

Proof of Fact 5: The family (f £ ) £ converges weakly to / in L^L^) — *, which with the strong 
convergence (Fact 3) of (u E * (f e ) e in L 2 X c -> Lj x , ensures the weak Lj x ^ convergence. Similarly, 

( TO o(/e7e))e converges weakly to mo/ in L^°(L^/ 3 ) — *, and Fact 4 ensures for instance that (u E ) E 
converges strongly in (L x ^ 2 ). 

Fact 6 [(1 + p £ ) (u E -kipe) ® u E ] converges weakly in hj x to (1 + p) u ® u. 

Proof of Fact 6: Once again, it's a " weak x strong " type of convergence here. First notice 

2 _ e i-e 

11 ~ 6 + 2 ' 

where 6* := 21/22. Hence, by Holder's inequality, we get the following interpolation result: L^°(L X ) (1 

2fi 2/0/11/2 2/0/11/2 

L t (L x ) c ^ Lj (L x ' ). This implies that both (u E ) E and (u E * ip £ ) £ are bounded in L t ' (L x ' ), and 
the product [(u £ -kip £ )®u £ ] £ (which converges almost everywhere to u<E>u, see Fact 4) is hence bounded 
in L^ 9 (L X 1 / 4 ). Since 9 < 1 and 11/4 > 5/2, we get the strong convergence of this product to u <g> u, 
in Lj (Ljj/ 2 ), and we already know that (p s ) e converges weakly to p in L^°(L X / 3 ) — 

Using all the weak convergences above, we end up with a local in time solution of our system, which 
concludes the Proof of Proposition [J|Tj □. 

Next subsection is devoted to the prolongation of the local solutions defined above in solutions defined 
on [0,T]. 

3.4 Energy estimate and global existence 

In order to prove the existence of global solutions to our system, a standard strategy consists in reproducing 
the previous step at time i* — e, obtain another local solution and paste it with the previous one, and so on 
and so forth. 

For such a strategy to succeed, one must ensure that the sequence of local times of existence does not 
decrease too quickly. But, as noticed in Rcmark l3.1( the local time of existence is bounded below by l/P™, 
which is (by definition) a non-increasing function of \\y/l + pi n Ui„||2 + ||/m||oo + ^hfin- We already know 
(by maximum principle) that ||/||oo < e 2t */i n , and this bound ensures that all possible extension will always 
satisfy ||/||oo < e 2T /i n , since the sum of all local times of existence does not exceed T (well if it does, we're 
done !). It is hence sufficient to prove that, for almost all t G J; n , 

M 2 f(t) + || Vl + p(*) u(t)|| 2 < M 2 /in + || v / l + p^u ln || 2 . (55) 

Indeed, such an estimate would propagate for each local solution and we may hence bound from below all 
the corresponding times of existence, which means that our strategy would end in a finite number of steps. 
In fact ([55)1 is a straightforward consequence of the following energy estimate: 
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Proposition 3.2. The solution built in Proposition ],']. 1\ satisfies for almost all t G J ln , 

\{M 2 f{t) + \\^irW)v{t)\\l H Ti)}+ / HVxu( S )||2 2( d s +| f f |u-£| 2 /d£dxd S 

1 y > JO ^ JO JPxl 3 

< ^{M 2 / ln +||Vl + PinU in ||^ (T3) }. (56) 

Proof of Proposition [3721 

We first use inequality for solutions of the regularized system, and add the integral in time of (|4"5|) x i 
to get 



\ [M 2 fS) + \\VT+W)Mt)\\l^)} + 1 IIVxu e ( s )||£ 2 



( T 3)ds 



3 '* 



K-£| 2 / £ d£dxd S 



JT 3 



< \ {vM 2 / in + \yj\ + p in UtaHjafpij} +Re{t), 

where 



Re{t) = l I I /e|u £ | 2 (l-7 £ (l))d|dxd S + 2 f f /^.U e ( 7e (^)-l)d|dxd S 

L JO JT 3 xt 3 Jo JT 3 xK 3 



Ri(t) 



/ / fs£. ■ (u e *<£ E - u £ )d£dxds 

JO JT 3 xR 3 



Using Lemma 13. 11 the maximum principle, and bound (|49|) . we notice that for a > small enough, (m a f e ) e 
is bounded in L^°(L^/ 2 ). Now recall that j E is chosen so that |1 — 7 E |(£) < l|£|>i/ E , hence 

which goes to with e since (u e ) e is bounded in L 2 (H X ) L 2 (L X ). This shows that (Rl) E — >0 uniformly 

£— >0 

on J; n , and a similar Proof applies for (i? 2 ) e . Then, we have 



R 3 e (t) 



[ [ (/ E £-v E )d£dxds, 

JO JT 3 xl 3 



where (v E ) e is bounded in L 2 (L X ) and converges to in L 2 X <-i L tx . For values of £ satisfying |£| < 
(||v E ||i + e)~ 1 / 2 , we simply use the maximum principle for (/ E ) E to see that their contribution goes to with 
e (uniformly in time), and for the large values of £ satisfying the opposite inequality, we handle them as we 
did for R\ and i? 2 , using the L 2 (L X ) bound. 

At this stage, we proved that (R e ) e converges to uniformly in time. Let us treat the other terms of 
inequality ([5d]) . We first use the classical estimates of weak convergence to get 



Avxu( s ; 

Jo 



lL2 (T 3)ds < Jim f HVxUellj^n^ds, 

e^O Jo 

M 2 f(t) < lim M 2 f e (t). (57) 

£->0 



Furthermore, adapting the proof of estimate (f5"4")) used in Fact 2, we obtain for almost all times t G J u 

\Wl + p e {t) u £ (t)|| 2 2(T3) ^ || v/l + p(t)u(i)|| 2 2(T3) , 
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so that it only remains to prove 

(I |u - £| 2 /d£dxcLs < Um f f |u e - £ | 2 / e d|dxds, (58) 

J JT 3 xR 3 e^oJo JT 3 xR 3 

in order to conclude the Proof of Proposition ^. 21 We first write 

/ f |u £ - £| 2 / £ d£dxds = f f |u £ | 2 /ed|dxd S + / M 2 f e (s)ds-2 f f £ ■ u £ / E d£dxds. 

Jo ii 3 xi 3 Jo Jt 3 xr 3 Jo Jo Jt 3 xR 3 

(59) 

Since (u e ) e converges strongly in L 2 X and (/ £ ) e converges weakly in L^ ^ — *, we have by Fatou's Lemma 
/ / |u| 2 /l|£|<„d£dxd S < lim / / |u e | 2 / £ d|dxd S , 

JO JT 3 xR 3 e^OJo JT 3 xM 3 

hence by monotone convergence, 

/ / |u| 2 /d£dxds < Um f f |u e | 2 / e d£ dxds. (60) 

JO JT 3 xR 3 £->OJo JT 3 xR 3 

Using another time Fatou's Lemma with estimate (|57[) shows that 



M 2 /(s)ds < lim / M 2 f e (s)ds. (61) 

e^O Jo 

On the other hand, since v e := u £ — u is bounded in L 2 (L X ) and converges to in L] x , we have as before 
(see the study of R^) 

f f /^.(iie-ujd^dxda— £0, 

JO JT 3 xR 3 

and handling large velocities as we did for R\ and i? 2 , we can show that 

/ (/£-u)l| € | >n d£dxds + sup / / (/ £ £ • u)l| € | >n d£dxds — > 0. 

JT 3 xR 3 £>OJo JT 3 xR 3 n->+oo 

When n s N is fixed, we have by weak LJ*^ ^ — * convergence of (f e ) e , 

[ ( (/ E £.u)l| € |< n d£dxds —»■ // (/ £ £-u)d£dxds, 

JO JT 3 xR 3 "^ + °°J0 JT 3 xR 3 

so that the three last convergences imply together 

// (/ e €-Ue)d£dxdfl — > // (/£-u)d£dxd S . 

JO JT 3 xR 3 »^ + °°Jo JT 3 xR 3 



(62) 



Using 0, (ED and ([121) in flSSJ, we get (JHHD and this ends the Proof of Proposition O □ 

This also concludes the Proof of Theorem 11.11 thanks to the strategy of prolongation of local solutions 
explained at the beginning of the Subsection. □ 



4 Appendix 

We present in the Appendix a few auxiliary Lemmas which are used in various parts of this work. We start 
with a standard variant of GronwalPs Lemma, that we recall for the sake of completeness. 
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4.1 A variant of Gronwall's Lemma 

Lemma 4.1. Let f £ < ^' 1 (R + ) fee a convex non- decreasing junction and a £ M. Consider (z, [0,4*[) 4/ie 
maximal solution of z' = f{z), z(0) = a on M + . Lei a : M + — > M + fee a continuous function such that 

a(t) < a+ ( f(a(s))ds. 
Jo 

Then, for all t £ [0, 4+[, 

a(4) < z(4). 

Proof of Lemma 14.11 

The usual (linear) Gronwall lemma shows that if u £ < rf 1 (M. + ) satisfies u(t) < / b(s)u(s)ds on [0,4* [, 

Jo 

then u is nonpositive on this interval. We have here by convexity (on [0,4*[) 

u(t) := a(t) - z(t) < [ [f(a(s))-f(z(s))]ds< f f'(a(s))u(s)ds, 



and since / is nondccreasing, u < on [0, 4*[. □ 

As a consequence, we get the 

Lemma 4.2. Let A, 7 > 0. TTie maximal solution of the Cauchy problem z' = Az 1+1 , z(0) = A is defined 
at least on [0, (7 A^ 1 )- 1 [. 

Proof of Lemma 14.21 

Indeed : 4 M> (A 1 — A 7 4) _ 1/7 7 is well defined for 4 < (7 J 4 7+1 )~ 1 , and solves the Cauchy problem. □ 
In next Subsection, we recall for the sake of completeness the following classical compactness Lemma: 

4.2 A compactness lemma 

Lemma 4.3. The injection x ^(L^) is compact. 



Proof of Lemma [ 

Let (f n )n be a bounded sequence of Hj x . 

Step 1: A = {t £ [0, T] : \\f n (t)\\ H i — > +00} is a null set with respect to the Lebesgue measure fi. Indeed, 
if it were not, then for any R > 0, we would have 

Rfi(A) < [ lim \\f n (t)f H i (4)d4. 

Fatou's lemma would then imply 

Rfi(A) < lim / ||/„(*)||hi < lim \\f n f H i , 

which is impossible since the right-hand side is finite. By definition of A, for all 4 in [0,T]\j4, the sequence 
(||/ r i(4)|j H i )n has a bounded subsequence. Since n(A) = 0, we can find a countable subset [0,T]\A that is 
dense in [0, T]. Let us denote this subset B = (4 p ) pg N- We can then extract a (diagonal) subsequence of 
(fn)n (sill denoted (f n )n) such that for all p £ N, the sequence (/n(4 p )) nG N is bounded in H x . 

Step 2: Since (dtf n )n is bounded in Lj x , (/„)„ is uniformly equicontinuous w.r.t. 4, with values in 

T 2 



1G 



Step 3: Since by Rellich's Theorem, the injection H x <->• L x is compact, we can extract (again, diagonally) 
a subsequence (still denoted) (/„)„ such that for all p, the sequence {f n {t P ))n converges in L x to some element 
f(t P ) G L x . 

Siep ^: / is uniformly continuous on B, with values in L x . Indeed, for any t, s in B, 

||/(f) - < ll/(*) - + ||/»(*) - /n(s)|| L| + ||/„(s) - f(s)\\ h% . 

The central term of the right-hand side goes to with \t — s|, independently of n because of the uniform 
equicontinuity obtained in Step 2. When t and s are fixed elements of B, the two other terms of the right-hand 
side vanish when n — > oo because of the extraction of Step 3. 

Step 5: f admits a unique continuous extension to [0, T], which is uniformly continuous on this interval, 
with values in L x . 

Step 6: f is actually the limit of (/„)„ in ^°(L X ). Indeed, for e > 0, choose 8 a corresponding modulus 
of equicontinuity for / and of uniform equicontinuity for (/„)„. Then pick (B is dense in the compact set 
[0, T]) a finite number t-y, . . . , of elements of B such that [0, T] C Ufl-Ji; — 8, U + 8[. Eventually, for any 
a e [0, T], if ex e]*j -8,ti + 8[, then 

ll/W - /n(cr)|| L| < - /(ti)!^ + \\f(U) - f n (U)\y x + \\f n (U) - f n (cr)\\ hi 

<2e + \\f(t i )-f n (t i )\\ hl , 

so that ||/(c) — /n(c)|lL 2 — 3e for n large enough. This ends the proof of Lemma T4. 31 □ 

4.3 Weak convergence of a product 

We present here a result based on the method used in [TU] . 

Proposition 4.1. Let q 6 [l,oo] and p 6 [1,3[. Consider two families (a e ) e and (b e ) e respectively in 
Lj(W x ' p ) and \j\ (L x ), with s < p* := -^z- Assume the weak convergences 

(a e ) e — >■ a, in L^(W X ' P ) — w-k, 

Q>e)e ~^ Q ^' *' n (L x ) — W ★ . 

If (dtb e ) e is bounded in Lf (H x m ) for some ra€Z, then, up to a subsequence, we have the weak convergence 

(a e b £ ) £ — v ab in the sense of measures, 

i.e. with test functions in < £f x . 
Proof of Proposition [47T1 : 

Let us first notice that for all e > a e b e € L] x (Sobolev injection) and that the sequences {a £ ) e and 
(b s ) e are bounded in the spaces in which they converge weakly. 

Step 1. We have 

a (b* ip n ) — > a b, in L] x strong . 
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Step 2. Since (d t b e ) E is bounded in Lf' (H~ m ), (& e *^r,) e is bounded in Wl' q ' (W x ' p ') so that (thanks to Rcllich's 
Theorem), for all fixed 77, (b e *<fn)e admits a (strongly) converging subsequence in hj (L x ), the limit 
being necessarily b-kip^ (this is due to the uniqueness of the weak— 7k- limit). In fact, we can choose (but 
we don't write it explicitly) a common (diagonal) extraction for all 77 after discretization (77 := 1/k). 
Since (a £ ) £ converges weakly in Lj (L x ), we eventually get, for all fixed 77, 

(a e (b £ * (p n )) £ ^ a(b-k ip v ) in x weak. 

Step 3. We shall use the following "commutator Lemma", the Proof of which is rather close to the usual 
Friedrichs Lemma (which is a key element of [TU]). 

Lemma 4.4. Under the assumption^ of Proposition \4~T\ if (Vjj)?) is a sequence of even mollifiers, 
then the commutator (convolution in x only) 

S e , v := a £ (b £ * tp v ) - (a e b £ ) * ip v 

goes to in L* x as 77 — > 0, uniformly in e. 

Proof of Lemma 14.41 

First recall the following standard fact : since (a e ) e is bounded in L' (W x p ) and s < p*, the sequence 
(rha e — a £ ) £ tends to in L*(L X ) as h — > 0, uniformly in e. We now write the following equality for 
the commutator 



S e ,r,(t,x) = o e (i,x) -a e (t,x-y) b £ (t,x - y) <p v (y) dy, 
whence thanks to Fubini's Theorem, 

\\S E , v \\i < II^IIl^ls') K(y)l \\ T y a e - a ellL?(L;) dy, 

which yields the desired uniform convergence, and concludes the Proof of Lemma 14.41 □ 
Step 4- We have 

(a £ b £ ) -k ifr. — a £ b £ — > in the sense of measures, 

uniformly in e (with a fixed continuous test function). Indeed, if 9 £ x , since ip^ is even, we know 
that 

((a £ b £ ) * Lp n — a £ 6 e , 9) = (a £ b £ , 9 -k ip v — 9) , 

and the right-hand side tends to with the desired uniformity because (a £ b £ ) £ is bounded in L* x , and 
(9 -k ip n — 9)^ goes to in LJ*^ (9 is uniformly continuous). 

Step 5. Write 

ab — a £ b £ = ab — a (6 ★ tp^) 

+ a (b * (p n ) - a £ (b £ ★ tp v ) 
+ a £ (b e * <p v ) - (a £ b e ) * <p v 
+ (a £ b E ) *tp v — a e b £ . 

Fix 9 £ "^Px- I n the right-hand side, line number i £ {1,2,3,4} corresponds to the Step i proven 
previously. We choose first 77 to handle (uniformly in e) all the lines of the right-hand side, except 
the second one. Then, we choose the appropriate e to handle the second line, thanks to Step 2. This 
concludes the Proof of Proposition ^. II □ 



1 The assumption on is obviously useless here. 
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Proposition 4.2. In Proposition \4-l\ when q < oo , the same conclusion holds assuming only a bound in 
(H~' rl ) (instead ofXj\ (H~ m ) ) for (dtb E ) E if, in addition, we assume the convergence of (a s ) s to a in L[ (L£) 
weak—k, for some r > q. 

Proof of Proposition [4721 

The Proof is identical to the Proof of Proposition 14.11 except for the second step (the only one using 
the bound on (dtb e ) e ). For this step, we use a (diagonal) extraction such that, for all fixed i], (b £ * ip v ) e 
converges strongly (and almost everywhere) to b*Lp n (but only) in L| (L£ ). Since (b E ) E (and hence {b e -kip 1 - j ) e ) 
is bounded in Lf (L£ ), with q' > r' > 1, we see that (b £ * <Pn)e converges to b* <p n strongly in U t (L^ ), and 
the added assumption of weak-* convergence for (a e ) e allows hence to get 

(a e (b e * tp-n))e a(b-k ip v ) in L* x weak, 
which ends the Proof of Proposition ^. 21 □ 
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